#!/usr/bin/env python3

import xarray as xr

f1 = xr.open_dataset('./mz.p1.360x180x30.h0.nc')
f2 = xr.open_dataset('./mz.p2x1.360x180x30.h0.nc')

l = 17
for k in range(len(f1.lev)):
	for j in range(len(f1.lat)):
		for i in range(len(f1.lon)):
			diff = f1.div[l,k,j,i] - f2.div[l,k,j,i]
			if diff != 0:
				print(i + 1, j + 1, k + 1, float(diff))
